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Entanglement in the ground state of a many-body quantum system may arise when the local 
terms in the system Hamiltonian fail to commute with the interaction terms in the Hamiltonian. We 
quantify this phenomenon, demonstrating an analogy between ground-state entanglement and the 
phenomenon of frustration in spin systems. In particular, we prove that the amount of ground-state 
entanglement is bounded above by a measure of the extent to which interactions frustrate the local 
terms in the Hamiltonian. As a corollary, we show that the amount of ground-state entanglement 
is bounded above by a ratio between parameters characterizing the strength of interactions in the 
system, and the local energy scale. Finally, we prove a qualitatively similar result for other energy 
eigenstates of the system. 
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I. INTRODUCTION 

A central problem in physics is understanding the 
ground-state properties of a complex many-body Hamil- 
tonian, especially the ground-state correlations. As an 
outgrowth of that interest, there has recently been con- 
siderable work on understanding the non-classical cor- 
relations in the ground state, that is, the ground-state 
entanglement. Some recent work on this problem, with 
further references^^ includes llliadlllESIllllllIl 
m m m El Ea. This work has been motivated by 
the remarkable recent progress in using entanglement as 
a physical resource to accomplish feats such as quantum 
computation and quantum teleportation^. 

In this paper we connect the phenomenon of ground- 
state entanglement to a well-known idea in condensed 
matter physics, that of frustration, which we now briefly 
review. More detailed introductions may be found in [iSj . 
A typical example of a frustrated spin system is shown 
in Fig. ^ It consists of a triangular arrangement of 
three spin-i particles, each pair being coupled by a clas- 
sical anti-ferromagnetic coupling (-t-JcTzCr^, with positive 
coupling strength J). The anti- ferromagnetic coupling 
means that neighbours prefer to be anti-aligned in order 
to minimize their interaction energies. However, a little 
thought shows that it is impossible for all three spins to 
simultaneously be anti-aligned with each of their neigh- 
bours. It is therefore not possible to simultaneously min- 
imize all three interaction energies, and the system is 
said to be frustrated for this reason. The ground state of 
the Hamiltonian is a compromise between the minimum 
energy states of the interaction terms. 
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t Electronic address: [nielsen^^hysics. uq.edu.au] 

URL: www . qin f o . org/people/ni elsen| 
1 See ha hj for reviews and further references. 



FIG. 1: A system containing three spin-i particles, coupled 
by a classical anti-ferromagnetic coupling {-i-Ja^az, with pos- 
itive coupling strength J) favouring anti-alignment. There is 
no way all the competing couplings can be simultaneously 
satished; for this reason we say the system is frustrated. 



Let us consider an analogous example in which frus- 
tration arises not from the difficulty of choosing simulta- 
neously compatible spin configurations, but rather from 
choosing simultaneously compatible bases for Hilbcrt 
space. For example, consider a system of two spin-^ par- 
ticles with Hamiltonian H — —g {a]. + a'^'j — a^a"^, where 
the superscripts indicate which spin the operators act on, 
and a^,ay, and cr^ are the usual Pauli spin operators. 
The ground state of this system arises as the result of a 
competitive process between minimizing the contribution 
to the energy from the local Hamiltonian, —g {a\. -f cr^) , 
and from the interaction Hamiltonian, —a\a'\. Of course, 
because these two Hamiltonians do not have common 
eigenvectors, the actual ground state cannot possibly 
minimize both simultaneously, and must be a compro- 
mise between the respective ground states of the local 
and interaction Hamiltonians. 

This example suggests a connection between the 
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ground-state entanglement and a generalized concept of 
frustration. If the interaction term in the Hamiltonian 
were turned off, the system would sit in an unentangled 
state — the ground state of the local Hamiltonian. As the 
interaction term is turned on, it causes the local Hamilto- 
nian to become frustrated. As a result, the ground state 
sits in a basis which is a compromise between the unen- 
tangled basis of the local Hamiltonian, and the basis for 
the interaction Hamiltonian. Provided the interaction 
was chosen appropriately, the result will be an entan- 
gled ground state. Furthermore, it is clear that the more 
frustrated the local Hamiltonian is by the interaction, the 
greater the potential entanglement in the ground state. 

The main result of this paper is a bound that makes 
these intuitive ideas quantitatively precise. Our paper 
thus illustrates a general idea discussed in [l|, [Tj, 120, 
I21I l22j| , namely, that quantum information science pro- 
vides tools and perspectives for understanding the prop- 
erties of complex quantum systems, complementary to 
the existing tools of quantum many-body physics. 

We begin in Sec. ^ by reviewing some basic material 
on quantitative measures of entanglement. In Sec. lIIII we 
prove a general, non-perturbative bound on the ground- 
state entanglement, relating it to the extent to which the 
interaction Hamiltonian frustrates the local Hamiltonian. 
We call this the "entanglement-frustration" bound. The 
proof of the bound is conceptually and mathematically 
extremely simple. Its interest lies in illustrating quanti- 
tatively a connection between two apparently disparate 
physical phenomena, and in the consequences which fol- 
low from this connection, to be discussed in later sections. 

In Sec. IIVI we apply the entanglement-frustration 
bound to an illustrative example. Using this example, 
we determine necessary conditions for the bound to sat- 
urate the ground-state entanglement. It is then shown by 
construction that it is possible to come arbitrarily close 
to saturation for all possible values of the ground-state 
entanglement, and we conclude that the entanglement- 
frustration bound is thus the strongest possible bound of 
its type. 

Aside from its intuitive appeal and immediate rele- 
vance, the entanglement-frustration bound has an ele- 
gant corollary described in Sec. El Intuitively, it is clear 
that the ground-state entanglement of a Hamiltonian 
H ~ Hl + Hi is small if the size of the interaction Hi is 
small compared with some appropriate local energy scale 
associated with Hi^. Indeed, it is straightforward to use 
perturbation theory to demonstrate a bound along these 
lines, valid in the limit when Hi is a small perturba- 
tion. The entanglement-frustration bound allows us to 
prove a general non-perturbative bound quantifying this 
intuition. This corollary is proved in Sec. Sec. IVII 
generalizes these results so that they apply to arbitrary 
eigenstates of the Hamiltonian, not just the ground state. 
This is done using methods quite different from those 
used in Sec. using a variant on a powerful theorem 
from linear algebra known as the Davis-Kahan theorem. 

The results in Sees. IIIIIVI provide a compelling pic- 



ture of how ground-state entanglement arises as the re- 
sult of frustration between competing local and interac- 
tion terms in the system Hamiltonian. Sec lVII eeneralizes 
some of these results to apply to other energy eigenstates 
as well. The paper concludes in Sec. lVIII with a discussion 
of some possible extensions to this work. 



II. BACKGROUND ON ENTANGLEMENT 
MEASURES 

To make our ideas precise we must introduce a quan- 
titative measure of the amount of entanglement in the 
ground state of a quantum system. A major focus of 
research in quantum information science over the past 
few years has been developing such a theory of en- 
tanglement^, and several good candidate measures ex- 
ist. We shall use a measure of entanglement introduced 
in [25I l26l| . For an n-body quantum system in a state V' 
this entanglement measure is defined by^: 

E{iIj) = 1~ max | . . . (g) V-n) p. (1) 

That is, -£^("0) measures the maximal overlap ■0 has with 
a product state ipi® . . .®^n of the n bodies making up 
the system. 

What makes a good entanglement measure? 

[25I [2^ investigated the properties of E{i}}) and found 
that it has many properties that make it a good mea- 
sure of entanglement. These properties include the fact 
that: (a) -£(0) can only decrease, never increase, under 
local operations and classical communication, i.e., it is 
an entanglement monotone; and (b) E{il}) is zero if and 
only if ■0 is unentangled, and otherwise is positive. In 
addition, an interesting connection has been found p7| 
between i?('0) and the theory of quantum algorithms, 
with E{il)) being related to the probability of success of 
an algorithm whose initial state is equivalent to i/;, up to 
a local unitary transformation. 



III. THE ENTANGLEMENT-FRUSTRATION 
BOUND 

The general scenario we consider is an n-body quantum 
system with Hamiltonian H = H^ + Hi. Hl is a local 
Hamiltonian consisting of single-body or local terms, and 
therefore has an eigenbasis of unentangled states. Hi 
contains all the remaining terms in the Hamiltonian, and 
is called the interaction Hamiltonian. 



^ See, e.g., l23l |24| for an introduction and further references on 
the theory of entanglement. 

^ N ote tha t this measure is a slightly rescaled version of that 
in |2,'ll2fill . but has essentially the same properties. In the present 
context the rescaled definition turns out to be easier to work with. 
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We let Eq be the global ground-state energy, i.e., the 
ground state energy of H, with \Eq) any corresponding 
ground state. Similarly Eq and Eq are defined to be 
the local and interaction ground-state energies, respec- 
tively, for Hl and Hj. We define the frustration energy 
of the system as Ef = Eq — Eq — Eq. The frustration en- 
ergy thus measures the extent to which the global ground 
state fails to simultaneously minimize the local and inter- 
action energies. It is easily shown from matrix eigenvalue 
inequalities that Eq > Eq + Eq, so Ef is always a non- 
negative quantity, and is equal to zero if and only if Hj^ 
and Hi have a common ground state. 

Our aim is to relate the amount of entanglement in 
the ground state, E{\Eo)), to the frustration energy, Ef. 
Of course, to relate the dimensionless quantity E{\Eo)) 
to Ef, which has units of energy, we require another en- 
ergy scale in the system. The relevant energy scale turns 
out to be associated with local excitations of the sys- 
tem. Suppose we decompose as Hi + H2 + . . . -I- iJ„, 
where Hj is the contribution to the local Hamiltonian 
from the jth system. Let AEj be the gap between the 
ground and first excited energies* of Hj. Let Ailment be 
the second smallest of these energies. That is, suppose 
we choose jQ,ji, . . . such that AEj^ < AEj-^ < • ■ •• Then 

AEcnt — ^Ej^ . 

Physically, AE'ont is the energy we need to put into a 
system with Hamiltonian H^ in order to cause an excita- 
tion from the ground state into an excited state of either 
system jq or system ji . It is thus the minimal amount of 
energy that we would need to put into the system in order 
to cause entanglement in the ground state, since merely 
exciting one system, while leaving the others alone, leaves 
the system still in a product state. 

Our result relating the ground-state entanglement to 
the frustration energy and Aii^cnt is the inequality: 



E{\E^)) < 



Ef 



AEr, 



(2) 



We call this the entanglement-frustration bound. This 
bound tells us that when the frustration energy is small 
compared with Ai^ont, there can't possibly be much en- 
tanglement in the ground state of the system. Thus, it 
is only systems in which the interaction and local terms 
substantially frustrate one another that it is possible to 
have a highly entangled ground state. 

The first step in the proof of the entanglement- 
frustration bound, Eq. is to prove that 



{Eq\Hl\Eq)-E^ <Ef. 



(3) 



Physically, this is just the obvious statement that the 
extent to which the local Hamiltonian is frustrated is no 



* Note for clarity: The term "gap" may be used in two different 
senses. Sometimes it means the energy difference between the 
ground state and the first excited state with a strictly higher 
energy. We use the term in the other sense, where the gap AEj 
is zero if Hj has a degenerate ground state. 



larger than the total frustration in the system. The proof 
is simply to split the frustration energy into a sum of 
contributions from the local and interaction frustration 
energies: 

Ef = {Eq\H\Eq) - E^ - El (4) 
= {{Eq\Hl\Eq) - E^) + {{Eq\Hi\Eq) - El) . 

(5) 

The inequality of Eq. Q now follows from the observa- 
tion that {Eq\Hi\Eq) > E^. 

The second step in the proof of the entanglement- 
frustration bound is to expand \Eq) in terms of the eigen- 
states lEj') of Hl, \Eq) = Ej ^il-^'j')- We assume that 
the local energies are ordered so that Eq < E^ < . . .. 
We now split the expansion of \Eq) into terms with ener- 
gies below Eq + Ai?ont, and into terms with energies at 
least Eq + ASont, that is. 



3=0 



(6) 



where (a) k is the largest integer such that Ej^ < Eq + 
AE^nt, and thus Ej^^-^ = + AE^nu (b) |£;_l) is a 
normalized state containing all the terms of energy at 
least Eq + AEcnt, and thus is orthogonal to the lower 
energy terms; and (c) 7 is the amplitude for \E±), and 
thus satisfies I7P = 1 — J2'j=o l^jl"^- 

For later use it is important to note that X]^=o ^j\Ej") 
is a product state, as all the terms \Ej) involve excita- 
tions of the same subsystem^ . Furthermore its overlap 
squared with \Eq) is given by J2'j=o 

Returning to the main line of the proof, from Eq. © 
we have 

k 

{Eq\Hi^\Eq) = ^ \ajfE^ + \jf{E^\Hi^\E^). (7) 
But E^ > E^, {E^\Hl\E^) > E^ + A£;e„t, and = 
fe 

{Eq\Hl\Eo) > ^|«,fi?o^ 



I {E^ + AE, 



Rearrangement of this inequality gives 



{Eq\Hl\Eq) - El^ > 



3=0 



AE^ 



(8) 



(9) 



System jo, to return to the notation used earlier in defining 
A Bent- 
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Combining Eqs. ^ and © we have 

Our desired result, Eq. will follow if we can establish 

that £;(|i;o)) < (l -Ej=o This follows imme- 

diately from the definition of the entanglement measure, 
Eq. 1^, and the observation we made earlier in the proof, 
that \Eo) and the product state J2 'j=o have over- 

lap squared X]j=o I'^iP- 

IV. APPLICATION AND SATURATION OF 
THE ENTANGLEMENT-FRUSTRATION BOUND 



Entanglement-frustration bounce 
Ground-state entanglement 
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FIG. 2: The ground-state entanglement and entanglement- 
frustration bound for the transverse Ising Hamiltonian IIH 
plotted against the parameter g. 



In this section we consider two separate but related 
issues. First, in Sec. IIV Al we apply the entanglement- 
frustration bound to an illustrative and physically rele- 
vant Hamiltonian, the two-spin transverse Ising model. 
This example is used to develop insight into the ques- 
tion of when the entanglement-frustration bound is sat- 
urated. Building on these insights, we analyse this ques- 



tion in more generality in Sec. IIVBI showing that the 
entanglement-frustration bound can be saturated for all 
possible values of the ground-state entanglement. Thus, 
there is a sense in which the entanglement-frustration 
bound is the best possible bound of its type. 



A. The two-spin transverse Ising model 

As an illustrative example, consider a system of two 
spin-i particles evolving under a transverse Ising Hamil- 
tonian 



H 



(11) 



In this model, the two particles are coupled magnetically 
along their z axes, and interact with an external magnetic 
field of strength g directed along the x axis. For the 
purposes of this example we take 5 > 0. The g < 
analysis is similar, but it simplifies the discussion to pick 
a definite value for the sign of g. 

Note that while the two-spin transverse Ising model is 
mathematically rather trivial, it has genuine physical in- 
terest. Furthermore, we will find that it is surprisingly in- 
formative as a way of understanding the conditions under 
which the entanglement-frustration bound is saturated. 
For these reasons we describe the results in some detail. 

Physically, 5 ^ is the strong coupling limit, where 
we expect the ground state to become quite entangled. 
We will see in detail below that it becomes maximally 
entangled in this limit, i.e., i?(|£'o)) \-, for our en- 
tanglement measure. In contrast, 5 ^ cx) is the weak 
coupling limit, and we expect the ground state should be 
a product state in that limit, E{\Eq]^ — > 0. 



The ground-state energy of Hll|) is easily found to be 
-£'0 — — -|- 4^2, and the ground-state is 

|i?o) = ^ ((25 + v/T+v) I + +> + !-->) (12) 

where = 1 + {2g -\- yT~M^)^ is a normalization con- 
stant, and |±) = (|0) ± |1))/V2. Note that |£:o> is 
in its Schmidt form, with largest Schmidt coefficient^ 

Ao - (2.9 + Vl + V) /V^- The ground -state entan- 
glement is given by 1 — Ag, which simplifies to 



£;(|£;o)) = - — ^ 



(13) 



To calculate the entanglement-frustration bound we 
must first split the Hamiltonian into a local and inter- 
action part, i/i = —gio'i + crl), and Hj = —a],a1. With 
these choices we find that — —2g and Eq = — 1. 
The two spin systems each have the same local energy 
spectrum with the gap between the ground and excited 
states being 2g, so we have ASont — 2(?. This gives the 
entanglement-frustration bound 



Ef 



AE, 



1 + 2g - y/l + 4g2 
2.9 



(14) 



A comparison of the quantities appearing in Eqs. 
and (|14|l is shown in Fig. [21 Both the ground-state en- 
tanglement and the entanglement-frustration bound de- 
crease sharply as g increases from 0. For these small 
values of g the bound is approximately double the entan- 
glement. As g increases further the ground-state entan- 
glement decreases rapidly to 0, while the bound decreases 



^ By contrast, if g < the largest Schmidt coefficient is Aq = 
1/VlV. This is the main difference between the g < and g > 
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to more slowly. The entanglement-frustration bound is 
clearly not very tight in this case, although the qualita- 
tive behaviour of the bound and the actual ground-state 
entanglement is similar. 

We can identify two reasons for the failure to sat- 
urate the entanglement-frustration bound in this ex- 
ample. First, in the language of Sec. IIIII the quan- 
tity {E±^\Hl\Ea) is strictly larger than Eq + Ai?ont- 

We see from Eq. that [E'j^) = | ), and thus 

{Ex^\Hl\E^) ^ E^ + 2AE^nt- It follows that 

{Eo\Hl\Eo) = XlE^ + (1 - XDiE^ + 2ASent) (15) 

and upon substitution into Eq. jSjl this gives 

Ef = 2A^e„t(l - A^) + (£;o|i^/|£;o> - 4 (16) 



excited state of Hl corresponding to some eigenvalue E^ 
then 



or 



Ef 



= 2E{\Eo)) 



{Eo\Hi\Eo) - Ej 

/S.Ec.nt 



(17) 



The entanglement-frustration bound is therefore at least 
twice the ground-state entanglement with this choice of 
for all values of g. 
The second contribution to the excess is the term 
{{Eo\Hi\Eo)-El)//\Eo^f Physically, this is the ratio of 
the frustration of the interaction energy to the local en- 
ergy scale. The excess sharply increases from for small 
g, and decreases slowly as g ^ oo. For g greater than 
about 2 the ground-state entanglement is close to and 
the entanglement-frustration bound is composed almost 
entirely of this excess term. 



B. Saturation of the entanglement-frustration 
bound 



When, if ever, is the entanglement-frustration bound 
saturated? We will show in this section that for all possi- 
ble values of E{\Eo)) wc can find a Hamiltonian H whose 
ground state has that amount of entanglement, and sat- 
urates the entanglement-frustration bound as closely as 
desired. 

Interestingly, it turns out that it is not possible to ex- 
actly saturate the entanglement-frustration bound except 
in the extreme cases E{\Eq)) ~ and E{\Eq)) = 1. How- 
ever, as we show in this section, it is always possible to 
saturate the bound to as good an approximation as de- 
sired. 

To see that exact saturation is not possible, consider 
the necessary condition for saturation {Ei^\Hi\Eq) = Eq 
identified in the previous section. This condition implies 
that |£^o) is a ground state of Hj, and therefore also an 
eigenstate of Hl = H — Hj. Entanglement in an eigen- 
state of a local Hamiltonian is only possible if there is 
an associated degeneracy. If \Eo) is a ground state of 
Hl then we conclude that Ai?cnt — 0, the entanglement- 
frustration bound is undefined, and so saturation cer- 
tainly does not occur. On the other hand, if \Eq) is an 



{Eo\Hl\Eo) - Eq — E^ — Eq. 



(18) 



But since \Eq) is entangled, by assumption, we must 
have Ej" — Eq > AE^nt- Combining this with the re- 
suh Ef > {Eq\Hl\Eq) - E^ gives Ef/AE,nt > 1. In 
contrast the maximum values of E{\Eq)) for qubits is ^, 
and more generally for pairs of d-dimensional systems 
it is 1 — i. We conclude that it is not possible for the 
entanglement-frustration bound to exactly saturate, ex- 
cept when E{\Eq)) = or 1. 

The above analysis, however, says nothing for \Eo) ar- 
bitrarily close to a ground-state of Hj , and in these cases 
it is possible that the bound approaches saturation. 

Before dealing directly with the issue of saturation, it 
is helpful to address another issue, the question of a how 
a given many-body Hamiltonian H is to be split into 
local and interaction parts. Consider, for example, the 
transverse Ising Hamiltonian H = —g{<jl. -\- cr^) — o\a'\. 
In our earlier analysis we set ~ —gial. + cr^) and 

However, there is a certain arbitrariness in the splitting 
into local and interaction Hamiltonians. From a math- 
ematical point of view, there is nothing to stop us from 
splitting H up as H = H'j^+H'j, where H'^ is any desired 
local Hamiltonian, and we simply choose H'j = H — H'j^ . 
So, for example, we could choose H'j^ = —gel, and 
Hj = —ga'^ — ala'^. The reason for this ambiguity is that 
while the class of local Hamiltonians is perfectly well- 
defined, there is no similar definition of what it means for 
a Hamiltonian to be an interaction Hamiltonian. Failing 
to have such a definition, we are free to choose how- 
ever we like, compensating by choosing an appropriate 
interaction Hamiltonian. 

This freedom to choose a splitting into local and inter- 
action parts is reflected in the fact that the entanglement- 
frustration bound holds for any choice of splitting H — 
Hl + Hj. Of course, while E{\Eq)) is not affected by the 
splitting chosen, the quantities AEcnt and Ef are. As 
a result the exact value of the entanglement-frustration 
bound depends on the particular splitting chosen. We 
will use this freedom in choosing a splitting to engineer 
saturation in the entanglement-frustration bound. 

Physically, of course, there is often a reason to favour 
one splitting into local and interaction parts over another. 
For example, if we regard the transverse Ising Hamilto- 
nian as a model of two magnetically coupled spins placed 
in an external magnetic field, then there is a clearly- 
defined physical sense in which —g{al. -t- cr^) ought to 
be regarded as the local term in the Hamiltonian and 
— (T^cr^ as the interaction term. 

However, the same model Hamiltonian may describe 
many quite different physical systems, and it is not at all 
clear that the splitting into local and interaction Hamil- 
tonians will necessarily be the same for all these physical 
systems. A priori it does not seem that the mathematics 
of quantum mechanics distinguishes any special subclass 
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Let us generalize this example further. Suppose H 
is an arbitrary bipartite Hamiltonian acting on two d- 
dimensional systems, with ground state Schmidt decom- 
position 



d-l 



\Eq) = X\aobo) + ^Xj\ajbj 



(21) 



where we have chosen labels so that A is the largest 
Schmidt coefficient. In order to ensure that Eq. (|19|l is 
saturated we choose a splitting of H with Hl as follows: 



0.25 
E(|E„>) 



FIG. 3: Comparison of the ground-state entanglement and 
entanglement-frustration bounds for two choices of splitting 
in the transverse Ising model. The solid line denotes the ideal 
case of saturation. 



of interaction Hamiltonians, and this makes it impossi- 
ble to define a unique splitting of H into local and in- 
teraction parties on purely mathematical grounds. More 
importantly, from our point of view, the entanglement- 
frustration bound holds for any splitting whatsoever, re- 
gardless of its physical (or unphysical) origin, and it is 
interesting to address the question of which splitting gives 
the best value for the entanglement-frustration bound. 

Let us return now to the question of saturation, and 
to a closer investigation of the example of the transverse 
Ising model considered in the previous section. In this ex- 
ample the decomposition of \Eq) into eigenstates of the 
local Hamiltonian is equivalent to the Schmidt de- 
composition, and the largest Schmidt coefficient is given 
by |(i?^|£'o)|- Furthermore the inequality 



{E^\Hl\Ej_) >E^ + AE, 



(19) 



is strict because \E±) = | ) is an excitation of both 

subsystems, whereas Eq -t- AE'ent is the energy of a single 
excited subsystem. The excess is therefore the energy gap 
of the remaining subsystem. 

On the other hand if we take advantage of the possi- 
bility of different splittings of H to choose — — ffcr^ 
then there is zero energy associated with an excitation of 
the second subsystem and Eq. ((T^ becomes an equality. 
The interaction Hamiltonian is determined by the choice 
of local Hamiltonian, Hj = H — Hl — —gc^ ~ '^I'^z^ ^^'^ 
we calculate a second entanglement-frustration bound 



E 



f 



AE,, 



yi + 4ff2 - ^1 + g2 

2.9 



(20) 



The two bounds Eq. H14() and Eq. H20I) are plotted against 
the ground-state entanglement in Fig. |21 It is clear that 
this second choice for H]^ provides a substantially tighter 
bound, as we expect. 



Hl = -7|ao)(ao| ® / 



(22) 



where 7 > is a parameter that will be chosen later in 
order to best saturate the bound. It is clear that Eq. H21() 
is an expansion of \Eq) in an energy eigenbasis of Hl, of 
the same form as used in Eq. |S1 and thus that 



{E^\Hl\E^ 







AE, 



(23) 



It follows that for this choice of local Hamiltonian, 



^(|£^o» 



Ef 



AE, 



{{Eo\Hi\Eo)-E^) /AE.nt, (24) 



i.e., the amount by which the entanglement exceeds the 
entanglement-frustration bound is composed entirely of 
the second term identified earlier in Eq . H17|l . 

To minimize this excess we choose 7 small and positive. 
Observing that Hj = H — we may do perturbation 
theory in 7 to show: 



Eq-{Eq\Hl\Eo)+0{-i^) (25) 
= (So|if/|So) +0(7'). (26) 



where we used Hi = H — Hl in the second line. Using 
this fact and the observation Ai?ont = 7, we have 



{Eo\Hi\Eo) - E^ 



AE, 



= 0(7). 



(27) 



Taking the limit as 7 — > we see that the entanglement- 
frustration bound approaches the ground-state entangle- 
ment. 

In summary, we have shown: 

Proposition 1. Let H be an arbitrary bipartite Hamil- 
tonian. Then there exists a local Hamiltonian Hl and 
corresponding interaction Hamiltonian Hi such that the 
entanglement-frustration bound derived from the splitting 
H — Hl + Hj is arbitrarily close to the ground-state en- 
tanglement of H . 

This shown that, in principle, the entanglement- 
frustration bound may be arbitrarily close to saturation 
for all possible values of the ground-state entanglement 
E{\Eq)). We therefore conclude that the entanglement- 
frustration bound cannot be strengthened without using 
more detailed knowledge of the system properties. 
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Our results show that saturation of the entanglement- 
frustration bound is always possible with an appropriate 
choice of splitting. They do not, of course, tell us what 
splitting ought to be used, except in the unusual situation 
where one knows virtually everything about the ground- 
state already, in which case one may as well calculate the 
ground-state entanglement directly. Thus the content of 
Propositionn]is not that we ought to expect to calculate 
ground-state entanglement exactly, merely by choosing 
the appropriate splitting for the Hamiltonian. Rather, 
Proposition n and the methods that lead to it, tell us 
that the entanglement-frustration bound is the best pos- 
sible, and provide some physical guidance as to how to 
choose the splitting into local and interaction Hamilto- 
nians in order to achieve the best possible values for the 
entanglement frustration bound. 



V. GROUND-STATE ENTANGLEMENT AND 
THE RATIO OF INTERACTION STRENGTH TO 
THE LOCAL ENERGY SCALE 

The inequality Eq. Q has a nice corollary that is easily 
proved, relating the ground-state entanglement to a ratio 
of the interaction strength with the local energy scale of 
the system. Suppose we define E^^^ to be the largest 
eigenvalue of Hj, and let E^^^ = E'^ax ~ be the to- 
tal energy scale for the interaction Hamiltonian, i.e., the 
difference between the largest and the smallest energies. 
It follows that 

Eo < {E^\H\E^) (28) 
= {E^\Hl\E^) + {E^\Hi\E^) (29) 
< E^+El,^. (30) 

Rearranging this inequality we obtain Ef < E^^^. Com- 
bining with Eq. ^ then gives 



E{\Eo)) < 



(31) 



The inequality Eq. (|31|l is an interesting result. Intu- 
ition, experience, and perturbation theory tell us that if 
we start with a local Hamiltonian and slowly turn on an 
interaction, the ground-state entanglement will depend 
on how strong the interaction is, compared with the lo- 
cal terms in the Hamiltonian, which tend to keep the 
ground state unentangled. Eq. H31() is a precise, com- 
pletely general statement of this intuition, a statement 
that holds even non-perturbatively. 



depend on how strong the interaction is compared with 
the strength of the local Hamiltonian. Of course, a sim- 
ilar intuition applies also for higher-energy eigenstates. 
Unfortunately, the strategy used to prove Eq. (|31|) cannot 
be applied directly to energy eigenstates other than the 
ground and most excited states^. The reason is that the 
proof of Eq. (|31|l relied on the entanglement-frustration 
bound, Eq. (|2Jl, and there is no natural analogue of this 
bound — or even a definition of frustration energy — for 
states other than the ground and most excited states. 

In this section we prove a bound validating this intu- 
ition for all energy eigenstates. The bound is proved in 
two steps. 

First, suppose A — B + C, where A and B are normal 
matrices. We will prove a general eigenspace perturba- 
tion theorem making precise the intuition that A and 
B have similar eigenspaces when C is sufficiently small. 
Our eigenspace perturbation theorem is a variant on a 
celebrated theorem of linear algebra, the Davis-Kahan 
theorem pj]^. 

A detailed discussion of how our eigenspace pertur- 
bation theorem compares to the Davis-Kahan theorem 
is given below. Summarizing, the major differences are 
that (a) our proof is simpler, (b) our conclusions are more 
powerful, but (c) our hypotheses are more specialized. 
For these reasons, we believe our eigenspace perturba- 
tion theorem is of substantial independent interest in its 
own right. 

The second step in the proof of the bound is to apply 
our eigenspace perturbation theorem to understand how 
the entanglement in an energy eigenstate depends on the 
relationship between the strength of the local and the 
interaction Hamiltonians. 

Let us begin with the eigenspace perturbation theorem. 

Theorem 1 (Eigenspace perturbation theorem). 

Let A, B and C he matrices such that A = B + C , with 
A and B normal matrices. Let a be an eigenvalue of A, 
and suppose Pa is any projector that projects onto some 
subspace of the corresponding eigenspace. (Pa may, for 
example, project onto the entire eigenspace.) Let (3 — {5} 
he some subset of the eigenvalues of B, and let Qb be pro- 
jectors projecting onto any subspaces of the corresponding 
eigenspaces of B. Define Q = X^he/jQb- Then 



(32) 



where S < T denotes a matrix inequality, i.e., T — S is 
a positive matrix, \S\ = VSS^, Aa = mint^fj \a — b\ is 



VI. HIGHER-ENERGY EIGENSTATES AND 
THE RATIO OF INTERACTION STRENGTH TO 
THE LOCAL ENERGY SCALE 

In Sec. we proved a bound, Eq. (|31|l . quantifying the 
intuition that when an interaction term is switched on in 
a many-body system, the ground-state entanglement will 



We only proved Eq. lijll for the ground state. An analogous 
result for the most excited state may be proved by applying 
Eq. I31i to the Hamiltonian —H. 
* For an account o f th e Davis-Kahan theorem, see Theorem VII. 3.1 
on page 211 of l2flll . and the surrounding discussion in Chap- 
ter VII of that work. 
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the distance from a to the set [3, and U is some unitary 
matrix. 

The interpretation of these inequahties in terms of 
eigenspace perturbation is perhaps not immediately 
clear. Rather than describe this interpretation imme- 
diately, we defer the description until after the proof of 
the theorem and a discussion of how this result relates to 
the Davis-Kahan theorem. 

Proof: We begin by proving the first inequality. Mul- 
tiplying A = B + C on the left by Pa and on the right by 
Qb, we obtain aPaQb = bPaQb + PaCQb, which may be 
rearranged to give 



PaQb 



PaCQb 

a — b 



(33) 



Observe that \PaQ\^ = PaQPa = EbPaQbQbPa- Substi- 
tuting Eq. (|33|l and its adjoint gives 



|pn|2 _ \ " PgCQbCPa 

PgCQbC^Pg 

A2 



b 



where we used |a — 6p > A^. Summing out b gives 



\PaQ\^ < 



\PaCQ\^' 
A2 



(34) 
(35) 

(36) 



The conclusion follows by using the operator monotonic- 
ity^ of the square root function, i.e., the fact that if 5 < T 
then VS <VT. 

To prove the second inequality in the statement of the 
theorem, it obviously suffices to prove \PaCQ\ < U\C\U^ . 
Note first that PaCQC^Pa < PaCC^Pa- But PaCC^Pa 
and C^PaC are positive operators with the same eigen- 
values, so there exists a unitary V such that PaCC^Pa = 
VC^PaCV^ < VC^CV'^. Putting these observations to- 
gether gives PaCQC'^Pa < VCtcyt^ from which it fol- 
lows that PACQC'iPa < UCC^W for some unitary U. 
The result now follows by using the operator monotonic- 
ity of the square root function. 

QED 

The conclusion of Theorem ^ has a nice implication 
in terms of matrix norms. Suppose ||| • ||| is a unitarily 
invariant matrix norm, i.e., |||C/S'V^||| — |||5||| for any uni- 
taries U and V. (Most of the familiar norms in common 
use in quantum information, including all the Ip norms, 
are easily shown to be unitarily invariant.) Using the po- 
lar decomposition we see that S = |S'|J7 for some unitary 
U, and thus Eq. (|32|l implies that 



\PaC 



^ \\\PaCQ\\\ ^ 
- Aa - 



for any unitarily invariant norm 1 1 1 • 1 1 1 . 

Let us compare the eigenspace perturbation theorem. 
Theorem n with the Davis-Kahan theorem. The Davis- 
Kahan theorem is as follows: 

Theorem 2 (Davis-Kahan theorem). Let A^B and 

C he matrices such that A ~ B + C , with A and B nor- 
mal matrices. Let a and [3 he subsets of the eigenvalues 
of A and B, respectively. Let P (resp. Q) project onto 
the space spanned by all the eigenspaces of A ( resp. B ) 
corresponding to elements of a (resp. (3). Suppose fur- 
thermore that a and (3 are separated by an annulus of 
width S in the complex plane, e.g., with a inside the an- 
nulus, and (3 outside the annulus. Then for any unitarily 
invariant norm III • |||, 



lll^'QIII< 



III^CQIII ^ IIICIII 



< 



(38) 



A, 



(37) 



There are three interesting differences between the 
Davis-Kahan theorem and Theorem ^ First, Theorem^ 
is more specialized than Davis-Kahan, in that it applies 
only for a single eigenvalue of A, not for multiple eigenval- 
ues. We have tried and failed to extend our proof to the 
more general case. A second difference is that TheoremQ] 
gives an operator inequality that implies the correspond- 
ing inequalities for unitarily invariant norms, but which 
is not implied by those inequalities. Finally, our proof of 
Theorem^seems to be substantially simpler than known 
proofs of the Davis-Kahan theorem. 

To better understand how Theorems Hand El relate to 
eigenspace perturbations, suppose that Pa projects onto 
a subspace Va spanned by a single eigenstate \a) of A, 
and Q projects onto a subspace Q spanned by eigenstates 
\b), b G (3. The norm |||Pa(3||| turns out to measure the 
orthogonality of these two subspaces. For example, in 
the special case when Q is a rank-1 projector, Q = \b){b\, 
we have 

|||P,Q||| = ||||a)(a|6)(6|||| = |(a|6)|||||a)(6|||| (39) 

which is proportional to the cosine of the angle between 
I a) and |5). (Note that ||| ||| is a constant inde- 
pendent of I a) and |6), due to unitary invariance of the 
norm.) Thus, Theorems ^ and |2] tell us that this cosine 
is very small (and thus \a) and \b) arc close to orthog- 
onal) whenever the ratio of the size of the perturbation 
|||C||| to the distance A^ is small. It follows that pro- 
vided |||C||| is sufficiently small, all the eigenvectors of A 
and B are nearly orthogonal, except for a single nearly 
parallel eigenvector. 

More generally, the singular values of PaQ are the 
cosines of what are known as the canonical angles be- 
tween the subspaces V and Q^". If |||PaQ||| is small then 



A rev iew of operator monotonicity may be found in Chapter V 
of I23. 



Fo r an introduction to the canonical angles, see Chapter VII 
of l2fll . especially the first section. We do not need to use any 
properties of the canonical angles in this paper. 
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the cosines of the canonical angles are small, and it can 
be shown that all vectors in V are very nearly orthogonal 
to all vectors in Q. 

Let us return now to the problem of bounding the en- 
tanglement in an arbitrary eigenstate \Ej) of a many- 
body Hamiltonian H. H is split into a local part, Hl, 
and an interaction part, Hj, as before. Our starting point 
is again the expansion of \Ej) in terms of the eigenstates 
\E^) of Hl. Associated to any local Hamiltonian we 
can identify some natural subspaces that contain no en- 
tanglement. These subspaces are spanned by a set of 
eigenstates \E^) related to each other by excitations or 
de- excitations of a single subsystem. Any superposition 
of such states factors into a product state, and for conve- 
nience we will refer to such a subspace as a product sub- 
space. Our use of this term should not be confused with 
the more general (and more common) use of the term 
product subspace, to mean any vector subspace contain- 
ing no entanglement; our use of the term is specific to a 
particular Hl, and refers to those subspaces spanned by 
sets of eigenstates \E^) which are all related by excita- 
tions or de-excitations of a single subsystem. 

We will see later that for each \Ej) there is a natural 
way to choose a corresponding product subspace from the 
eigenstates of Hl. For now let K. be any such product 
subspace and expand \Ej) in the energy eigenbasis of Hl 
as follows 



J2 ak\E^)+j\E, 



k, \Ef;)eK 



where |i?_L) is orthogonal to all states in /C. 
from Eq. (J) that 



(40) 



It follows 



E{\E,)) < 1 



fc 



I7l 



(41) 



Our strategy is to apply Theorem ^ to obtain a bound 

on |7p = l-EfeKP- 

Define Pj to be the projector onto \Ej). We're trying 
to bound the amplitude squared of the component 
of \Ej) orthogonal to /C, so let JC± denote the subspace 
spanned by all eigenstates ji?/") of Hl not in K, and de- 
fine Qk.^ to be the corresponding projector. Theorem ^ 
implies that 



\PjQk 



< 



\\\Hi\\\ 



= mm I 



\E, 



II is related 



El 
to 



where A^x 
show how 

It is easily seen from Eq. ()40|l that Qk. 
and so 



(42) 

Next we must 
the entanglement 

\E,) = \E^_) 



W^QkAW = ||||£;,)(£;,|Q^J|| (43) 

= \l\\\\\E,){EA\\\. (44) 

As remarked earlier, the value of ||| |w)(w| ||| for any nor- 
malized vectors \v) and \w) is a constant that depends 



only upon the norm ||| • |||. Without loss of generality 
we may assume that ||| |'y)(w| ||| = 1, since multiplying 
a unitarily invariant norm by a constant gives another 
unitarily invariant norm. We will say any norm satisfy- 
ing this condition is normalized. (Examples of normal- 
ized unitarily invariant norms include the operator norm 
||j4|| ~ ^'^\>\\\v)\\=\ the Hilbert-Schmidt norm 

PII2 = ^/MX4t).) 

Assuming that 1 1 1 • 1 1 1 is normalized we see that 



and it follows from Eqs. iglj and (021 that 



\Hr 



(45) 



(46) 



In Appendix^we show that for any normalized unitarily 
invariant norm ||| • ||| we have 115*11 < |||S'|||, where || • | 
is the operator norm and S any operator. The strongest 
form of the bound is therefore: 



E{\E,)) < 



(47) 



Different choices of the product subspace JC provide us 
with a different bound in Eq. 147|l . Ideally we would like 
to choose K. so that the quantity Aj-^x as large as 
possible. If Ej, or a good approximation to Ej, is known 
then we would ensure that /C contained |i?^) where \Ej^ — 
Ej I is minimal. More typically Ej is unknown, and this 
is not possible. However, there is still a natural way for 
us to choose /C. Importantly this choice also allows us to 
obtain a lower bound for AEj^jc^^ in terms of relatively 
simple quantities that depend only Hl and Hj, not on 
typically difficult-to-calculate quantities associated with 
the total Hamiltonian, H. 

Let {Ej") be the j-th excited eigenstate of the local 
Hamiltonian. We choose the product subspace JC so that 
the expression 



AE. 



j,cnt 



Ej 



Et\ 



is maximized. AiJ^ ont is a generalization of AE^ 



(48) 



Sec. mil in that it is the energy required to excite or de- 
excite at least two subsystems from the state \Ej"). Note 
that the calculation of AEj^cnt is tedious, but in prin- 
ciple straightforward provided that the energy spectrum 
of Hl is known: simply enumerate the possible product 
subspaces given the spectrum of Hl (a long, but finite 
list), and then calculate the minimum by inspection. 

Now for each \Ej^) G JC± we have by the triangle in- 
equality 



\Et 



E, 



> 



\Et 



> AE. 



'j.cnt 



-\Ef 
\E^ - E, 



E, 



(49) 
(50) 



Furthermore it is straightforward to show that \Ej 



< I-E^inaxl and so 



hX± = nun 



Ei 



Ei > AE. 



\Ei 



(51) 
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Substituting into Eq. (|47|l we obtain a result in terms of 
the spectrum of and the strength of Hj alone. 

Proposition 2. Let H ~ + Hj with H]^ a local 
Hamiltonian, and suppose A^^cnt > l-E-i^iaxI- Then the 
entanglement in the jth excited eigenstate \Ej) of H , as 
measured using the definition of Eq. is bounded above 
by: 

E{\E,)) < 2 (52) 

{AEj^ent - I-B^axl) 

Noting that I-E^jaxl — 11^/ 1 1 t^iis can be restated in a 
slightly weaker but perhaps more elegant form, supposing 
A,,c„t> 

Ei\E,))< M 1 (53) 

" - (Ai?,-e„t-||i//||)' (^ Ag.e„, 

Eqs. H52|) and H53|l confirm and quantify our intuition 
that when the non-entangled energy scale associated with 
\Ej) is large compared to the strength of the interaction 
Hamiltonian we expect little entanglement in the excited 
state \Ej) of the total Hamiltonian. 

Eq. H52(l should be compared with the earlier result 
Eq. H31|l for the ground-state entanglement. We see that 
the present result is equivalent to the earlier result, ex- 
cept for the presence of the term — jS'^axI i'^ the denomi- 
nator of Eq. (|52l) . which makes the present result weaker. 

VII. CONCLUSION 

We have introduced the frustration energy i?^ as a 
measure of the degree of frustration between local and 
interaction terms in the Hamiltonian H = + Hj of 
a many-body quantum system. This measure, when re- 
lated to a local energy scale, allowed us to derive the 
entanglement-frustration bound on the ground-state en- 
tanglement in the system. A novel feature of this bound 
is that it depends only on spectral properties of the 
Hamiltonians H,Hl and Hj. Ground-state entangle- 
ment properties can therefore be easily inferred directly 
from the spectra alone. 

The entanglement- frustration bound has, in turn, been 
used to prove a bound, Eq. (|31l) . relating the ground- 
state entanglement to a ratio of the strength of the in- 
teractions and an appropriate local energy scale. This 
bound involves only the eigenvalues of the local and in- 
teraction Hamiltonians, which are typically much easier 
to calculate than the eigenvalues of the full Hamiltonian, 
and thus this bound is more likely to be useful in prac- 
tice. A similar bound for an arbitrary energy eigenstate 
is proved in Eqs. and 

Ultimately it would be useful to have many powerful 
general techniques enabling us to infer ground-state en- 
tanglement properties of a Hamiltonian by considering 
the interplay between its constituent terms. This is not 



always easy. For example, consider the following system 
of three spin-i particles 

H = gaHA + gbHs + gcHc + Hab + Hbc (54) 

where A, B, C label the three particles. Ha, Hb, He are 
local Hamiltonians, Hab, Hbc are interaction Hamil- 
tonians on the appropriate subsystem, and gA,gB,gc 
control the respective strengths of the local Hamiltoni- 
ans. The bound Eq. H31|) derived from the entanglement- 
frustration bound tells us that if gi, is relatively large 
then there is little entanglement between particle B and 
the rest of the system AC. From this we may deduce 
that if there is any entanglement in the ground state 
then it must be between particles A and C. To some ex- 
tent, then, the entanglement-frustration bound allows us 
to determine the distribution of entanglement. In cases 
where all three local energy scales are small compared 
to the interactions, however, we are unable to directly 
deduce anything using the techniques in this paper. 

Throughout this paper we have defined frustration to 
occur when it is not possible to find a simultaneous 
ground state for some local and interaction part of a 
Hamiltonian. This is based on an analogy to the useful 
definition of frustration, which involves competition be- 
tween interactions, as discussed in the introduction, and 
illustrated Fig. ^ (An insightful review of classical and 
quantum frustration in this sense may be found in |3fl].) 
Both these points of view suggest interesting extensions 
of the investigations in the present paper. 

For example, we believe that quantum frustration sug- 
gests interesting parallels with the phenomenon of en- 
tanglement sharing 31] which places restrictions on the 
distribution of entanglement amongst many particles. In 
particular, we expect non-trivial distributions of entan- 
glement in the ground state of two overlapping interac- 
tions. For example, consider a Hamiltonian acting on 
three spin-i particles as before 

H = Hab + Hbc (55) 

and suppose that Hab smd Hbc have non-degenerate, 
maximally entangled ground states. It is impossible for 
entanglement to be distributed in a way that would pro- 
vide a ground state for H that is a simultaneous ground 
state of Hab and Hbc- The system is therefore nec- 
essarily frustrated. We might ask what happens to the 
ground-state entanglement distribution in systems such 
as this, and whether there are any properties of the con- 
stituent Hamiltonians that allow us to prove quantitative 
bounds relating the distribution of two-party, GHZ-type 
and W-type entanglement in this system. 
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APPENDIX A: UNITARILY INVARIANT 
NORMS AND THE SUP NORM 

In this Appendix we prove that for any normaHzed 
unitarily invariant norm ||| • |||, \\S\\ < \\\S\\\, where || • || is 
the usual matrix norm. This result is an easy corollary 
of the follo wing theorem, proved as Theorem 3.5.5 on 
page 204 of 0(: 

Theorem 3. Let \\\-\\\ be a unitarily invariant norm on 
the space of n x n matrices. Then there exists a compact 



set M in 3?" whose elements are decreasing sequences of 
positive real numbers, and such that 

n 

|||5||| = maxVm,a,(^), (Al) 

mGAf ^ — ^ 

where <yj{S) are the singular values of S, arranged into 
decreasing order. 

For a normalized unitarily invariant norm we have 
III |s)(s| III — 1, where |s) is any normalized state. It fol- 
lows from the theorem that there exists m g M such that 
TOi = 1. It follows that for any S, \\\S\\\ > ai{S) = ||5||, 
as we set out to prove. 
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